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Abstract. 

It is very common to ignore the non-local bulk effects in the study of brane- 
world cosmologies using the brane-world approach. However, we shall illustrate 
through the use of three different scenarios, that the non-local bulk-effect Vuu 
does indeed have significant impact on both the initial and future behaviour of 
brane-world cosmologies. 
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1. Introduction 

There is much interest in cosmologies derived from alternative theories of gravity. 
String theory, being one of these alternatives has motivated many different viable 
cosmologies. In particular there is a class of cosmologies motivated by string 
theory in which matter fields are restricted to a 3-dimensional brane embedded in 
1 + 3 + d dimensions (the bulk), while the gravitational field is free to propagate 
in the d extra dimensions [fl El El ED El El- I n fact, Randall and Sundrum [7| E] 
have been successful in recreating conventional Einstein gravity at low enough 
energies in their second brane-world proposal here-after labelled RS2. There 
has been many studies of cosmologies derived from these brane-world scenarios. For 
example, the so-called perfect fluid Friedmann brane-world models were investigated 
in EH El El E3 EI ED EH and more recently, people have investigated anisotropic 
brane-world models [El CEl 113 EHl HHI EDI EH E21 ESI EH E3 See [2E1 E3 EHI ESI for 
reviews of brane-world cosmology. 

There are essentially two approaches in studying brane-world cosmologies. Firstly, 
the bulk-based approach |51 1101 1X1*1114) in which one first solves the higher dimensional 
field equations and then attempts to determine whether a lower-dimensional brane- 
world exists and subsequently studies the cosmological dynamics on the brane. 
Secondly, the brane-world approach in which one solves the projected field equations 
on the brane |3(J1 131) and essentially ignores the extension of the solution to the bulk. 
Both of these fundamentally different approaches have their benefits as well as their 
drawbacks (See [T7] for additional arguments). 

One of the inherent problems of studying brane-world cosmologies using the 
brane-world approach is that observers on the brane cannot possibly determine all the 
effects resulting from the bulk. In particular there exist non-local effects (U, Q^^V^) 
coming from the bulk that affect the local dynamics of observers on the brane |31|. 
It is not at all well understood what the non-local effects from the bulk have on 
cosmological observations in the brane. In most papers that study RS2 brane-world 
cosmologies, one of the following assumptions is often made: 

• U = = V^u = 0, in which case the bulk is conformally flat and there are no 
non-local bulk effects affecting the dynamics on the brane. 

• U ^ but <2 M = V^u = 0, which are necessary conditions for the existence of 
a Friedmann Brane-world while the remaining non-zero component of the non- 
local effect on the brane, U, behaves as if it were a radiative fluid (assuming it is 
positive) , or can be interpreted holographically as the dual conformal field theory 
living on the brane |32l 133) . 

• U 7^ 0, Qfj, = and a^V^u — 0, a common constraint in the investigation 
of anisotropic but spatially homogeneous brane-worlds, chosen primarily for the 
ease in analyzing the models in a self consistent manner |19l 120) . 

There has been some additional work on determining the effects of a non-zero V^v on 
the behaviour of brane-world cosmological models |19l 1201 l2~2l 1251 134) . However, for 
the most part, all the assumptions above are essentially adhoc and not well founded. 
Below, we will make an alternative but reasonable assumption on the form of the 
unknown V^. 

A note with regards to notation. Tilde quantities denote quantities in the five 
dimensional bulk while non-tilde quantities are quantities on the four-dimensional 3- 
brane. The upper case latin indices range from to 4 while the lower case greek 
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indices range from to 3. Round (square) brackets enclosing indices represent (anti-) 
symmetrization. 



2. The Effective Field Equations on the Brane 

We shall investigate brane-world cosmologies using the brane-world approach. 
Shiromizu et al. jjQj have developed an elegant covariant approach to determine 
the effective four dimensional field equations induced on a single brane embedded 
in a five dimensional bulk using the Gauss-Codazzi equations, the so-called Israel 
junction conditions |35j and the usual Z2 symmetry. Maartens |31| further developed 
the ideas of Shiromizu et al. by extending earlier work of Ellis and MacCallum 
[361 137) by determining the dynamical equations describing braneworld cosmologies. 
The resulting equations are a modification of the standard Einstein Field equations, 
with new terms transferring bulk effects onto the brane. Assuming a five-dimensional 
bulk spacetime containing a 3-brane, the field equations are Einstein's equations with 
a (negative) bulk cosmological constant A and brane energy momentum Tab, as a 
source 



Gai 



-kgAB + 5(x){-^9ab + Tab} 



where x = yields the location of the 3-brane and k 2 — 8n/Alp. Let g M „ be the 
induced metric on the brane, then the effective field equations on the brane become 



where 



G fu At/^jy 
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k T^n 



« 2 = tt^> A = 6^ , A = ^ 
M p 2 K 4 Mp 3 
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(1) 



(2) 



The bulk corrections to the Einstein equations on the brane manifest themselves 
in two forms: firstly, any matter field on the brane contributes to local quadratic 
energy-momentum corrections via the tensor S^ v defined as 



O _ rri arp 



(3) 



Secondly, there are nonlocal effects from the free gravitational field in the bulk, 
transmitted through a projection of the bulk Weyl tensor onto the brane. If n A 
is the spacelike unit normal vector to the brane then we define 



CABCDn c n D 



which is both traceless and symmetric and has no components that are orthogonal to 
the brane. We then take 

lim £ A b = £av^J> B 

X^O An 

to define the non-local effect of the bulk gravitational field on the brane |3T]. Given 
a timelike congruence u^, on the brane, the bulk correction, can be decomposed 
into an equivalent non-local energy density on the brane, U, a nonlocal energy-flux on 
the brane, Q^, and a non-local anisotropic pressure on the brane, V^, via 

4 r 



1 



U( Ufl u u + -h^y) + TV + 2Q 



(4) 
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where 



4 



u = - = e^uv, (5) 



A 1 

K ' 

4 



Q, = (|) (6) 

?V = - (|) 4 ^<^>- (?) 

Where /i p „ = + u^u v projects orthogonal to on the brane and the angled 
brackets denote the projected, symmetric, and trace free part via 

v<»> = h;v„, w <flv> = [h { «h v f - h^h^w^ 

(See for complete details.) 

Due to the nature of the matter in the bulk (a cosmological constant) and the Zi 
symmetry, it can be shown that the brane energy momentum tensor on the brane is 
separately conserved, that is 

r*t„ = o. (8) 

It then follows from the contracted Bianchi identities on the brane that 

In general, equation © yields evolution equations for U and <2 M (see j3J) for 
further details). Unfortunately, there are no evolution equations for the non-local 
anisotropic pressure, "P M „. Therefore, generically the projection of the 5-dimensional 
field equations onto the brane does not lead to a closed system of equations. 



3. Governing Equations 

In order to illustrate the possible effects of V^u on the dynamics of a braneworld 
cosmological model, we first make some simplifying assumptions. 

• Firstly, it is common to assume through fine tuning of the brane tension A (a 
la Randall Sundrum) that the effective cosmological constant on the brane is 
zero. Here, for illustrative purposes only, we shall also assume that the effective 
cosmological constant on the brane is zero, that is, A = 0. 

• The simplest geometry for the braneworld, that permits a non-zero V^y is that of 
a Bianchi type I spacetime geometry. That is, we will assume a Bianchi I brane 
with line element of the form 

ds 2 = -dt 2 + [ ai {f)] 2 dx 2 + [a 2 {t)] 2 dy 2 + [a 3 (t)] 2 dz 2 . (10) 

• For simplicity, we shall assume that the matter on the brane is in the form of a 
perfect fluid with four velocity orthogonal to the surfaces of homogeneity (non- 
tilted). The energy- momentum tensor for a perfect fluid on the brane is given 
by 

= (WfiUu + phftu (11) 

where w M is the fluid 4- velocity and p and p are the energy density and isotropic 
pressure respectively. 
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It then follows that the Hubble expansion scalar, acceleration vector, and 
anisotropic shear tensor are of the form 

21 _l_ 21 _i_ 2± 
oi a 2 as. 



u tl - v u v = 

hi 



1 



-diag 0, 2 



0,2 a 3 



a 2 



2- 

a-3 «2 



Ol 

«1 



03 

a-3 ' o 3 



21 



«2 
«2 



ai a 2 

where an over-dot represents a derivative with respect to t. 

With the assumption of a perfect fluid source on the brane, the local bulk 
corrections, (see equation |3J) , to the effective four-dimensional field equations reduces 
to 



S,, 



1 



f 



j iiv - j^P u^u v + —p (p + 2p) h^y . 

Both the local and non-local bulk corrections mentioned above can be consolidated 
into an effective total energy density, pressure, heat flux and anisotropic stress as 
follows. The modified Einstein equations can be rewritten in the standard Einstein 
form with a redefined energy-momentum tensor: 



(12) 



n — k z t 

'-'III' — ru - 1 I 



2T-itotal 



where 



T, 



total 



i 



(13) 
(14) 



K 

K 2 ' 

is equivalent to an imperfect fluid energy-momentum tensor with the total equivalent 
energy density, pressure, heat flux and anisotropic pressure due to a perfect fluid 
source on the brane 



total 



total 



total 



total 
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K 4 


V 


K 6 


12 



-p (p + 2p) + \u 



6 ^ 



K 

—Zrpv- 
K 



(15) 

(16) 

(17) 
(18) 



The effective Einstein field equations (|13I14|) yield the following dynamical 
equations describing the evolution of the Bianchi I brane-world (Hll EHj : 

H = -JBP-^o*- — (p total + 3iO total ), 



3 6 



' //// 



,J1 ^total o tt2 Jl 
K p = 6H — (7 , 

«V° tal = 0, 



(19) 

(20) 
(21) 
(22) 



where p total , pta**^ gt° tal and 7r'° tal are defined in equations (|15I18|I and 
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From equations (assuming p — (7 — l)p) we obtain the following conservation 

equations 

p = -3-yHp, (23) 
U = -mU-a^V^, (24) 
hiP;. 6 = 0. (25) 
We define new dependent dynamical variables of the form 



V3^ 3tf 2 ' 

n o _ 3ff 2 f « a pY 

"~ K 2 A3ff 2 ' A "2A« 2 V3i? 2 y 



(26) 



and a new independent variable 

dt _ 1 
~&r ~ ~H' 

We also define a dimcnsionless variable V to represent the nondocal bulk corrections 
arising from the quantity V^v 

v 6 V3 0^ 

k2X sftH\[\o~o^ k2X aH2 



The resulting system of differential equations decouples into an evolution equation for 
H, (which we use to define the deceleration q, ) 

H' = -(q + l)H, (28) 

so that 

q = 2E 2 + n u + i(37 - 2)0 + ^(67 - 2)fi A , (29) 

and a reduced dynamical system describing the dynamics of the Bianchi I perfect fluid 
braneworld with nondocal bulk effects 

£' = T,(q-2)+V, (30) 

O' = Sl[2q - (3 7 -2)], (31) 

n' A =n A [2g-(67-2)], (32) 

0^ = n w (2g-2)-2SP, (33) 

where ' denotes differentiation with respect to t. Equation l|21|) yields a constraint of 
the form 

1 = s 2 + n + n x + n u , (34) 

that can be used to globally eliminate one of the dynamical variables. 

There is no evolution equation for the quantity V. This is not unexpected since 
observers on the brane are not expected to be able to determine all characteristics of 
the bulk (see OEHI)- In order to proceed further, we must make some assumption 
as to the form of V . In order to obtain a closed system of equations, it would not be 
unreasonable to expect that V is a differentiablc function of the remaining dynamical 
variables, that is, 

v = 7>(E,n,fi A ,n w ). (35) 
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which clearly satisfies the conservation equation 1)25(1. With this assumption, the 
value of q (see equation 12911 at the equilibrium points of the dynamical system (130(1 
is constant. Provided g ^ -1 at the equilibrium point, we are able to solve equation 
1(28(1 and obtain the solution H(t) = Hot -1 . With this, we are now able to find the 
solution of the remaining variables, p = pot~ 2 , — (a^Qt" 1 , U = Uot~ 2 , and 
V^v = (V^ jot^ 2 where H , p , (07^)0, Uq, {P^)o are all constants. It is worth noting, 
that the resulting braneworld spacetime represented by these equilibrium points is 
geometrically self-similar j^Hl I40j . (If q = —1 at an equilibrium point, then the 
resulting solution has H, p, o~,U, all equal to constants and consequently the resulting 
spacetime is not self-similar.) The property of asymptotic self-similarity is one that 
has been investigated by many in the context of general relativistic cosmologies, see 
(391 1401 l4~T] . Here we observe, that our assumption on the form of V is a necessary 
condition in order for the equilibrium points of the resulting reduced dynamical system 
to represent self-similar cosmological models. Now it is not at all clear what the actual 
form of this function should be, but we shall investigate a number of cases. 

• V = Vo, a constant: This case includes the scenarios — studied by 
0[T2[Tim[T3[Iinimi23EI3EIland V^a^ = [HOni as special subcases 
oiV = 0. 

• V = "P^X", a > 0: This ansatz guarantees that the necessary condition 
■p(S, ri, ^£y)|e=o = is satisfied for the existence of an isotropic equilibrium 
point. 

• V = V^fl^, (3 > 0: This ansatz is perhaps more desirable than the others in 
that both V and fly are directly related to the Weyl curvature in the bulk, while 
all the other dynamical variables are inherently on the brane. This assumption 
is similar to that used in [221 12] ■ 

With an assumption on the form of V, the evolution equations for the variables 
X — [£, f2, flu] yields a dynamical system of the form X' = F(X) suitable 
for a qualitative analysis using dynamical systems techniques |37| . In short, the 
procedure requires one to find equilibrium solutions of the dynamical system, that 
is, points X = X$ such that F(Xq) — 0. If the eigenvalues of the linearization of 
F(X) at X = Xq have non-zero real parts, then the nonlinear system X' = F(X) 
behaves qualitatively the same as the linear system X 1 = DF(Xq)(X — Xq) in a local 
neighborhood of X = Xq. If the real parts of the eigenvalues of DF(Xq) are all 
negative (positive), then the equilibrium point is said to be stable to the future (past) 
as all orbits in a local neighborhood of the equilibrium point have this equilibrium 
point as a future (past) asymptotic state, and the point is consequently called a 
sink (source) . If the eigenvalues have both positive and negative real parts, then the 
equilibrium point is called a saddle and is subsequently unstable to both the past and 
to the future. 

4. Qualitative Analysis 

4-1- V = Vq 

Equations The dynamical system (IrtOJ) is invariant under the transformation 
(■p,E) — * (— V,— £). This implies that we can restrict ourselves to the scenario Vq > 
without loss of generality. The required analysis becomes more manageable if we also 
redefine the parameter Vq — f(P) — P(P 2 — 1) where / : [1, 00) — > [0, 00) is a 



Brane-world Cosmologies with non-local bulk effects 



homeomorphism. Equation (|34[1 can be used to eliminate the variable fly from the 
dynamical system (|3Ull to yield the following three dimensional system 

£' = E(q- 2) + P(P 2 - 1), 

Cl' = n[2q- (37-2)], (36) 
^ = n A [2g-(6 7 -2)], 



where 



9=l + S 2 + i(3 7 -4)0+i(6 7 -4)fi A . (37) 
The equilibrium points and their local stability are discussed below. 

4-1-2. Equilibrium Points Bianchi I perfect fluid model (BI p /): 

3(2-7) ' 9(2-7) 2 (37-4)' 2 

This equilibrium point represents an anisotropic Bianchi I perfect fluid cosmology. 
Because 1 < P < 00, we observe that E > 0. This point is physical only when fl > 0, 
which implies that P and 7 must also satisfy 

8P 2 ( p2_ 1)2 

1 - 9 (2-7) 2 (37-4)^°- (38) 

For 7 < 4/3 the constraint is satisfied for all values of P. However, for 4/3 < 7 < 2 
the constraint is essentially an inverted cubic function of P 2 with up to three roots. 
It can be shown that P 2 — 1 is a local max for all values of 7, and that the largest 
real root is P 2 = Pi where 

1 



Pl 4 



-(3 7 - 8) + \/(3 7 - 8) 2 - 36( 7 - 2)< 



Therefore the constraint is satisfied if P 2 is less than the value of the largest root, 
i.e, P 2 < P\. In the limit as the parameter P — * 1 this equilibrium point becomes 
isotropic and represents a zero curvature perfect fluid Friedmann-Robertson- Walker 
cosmology. The exact solution at this equilibrium point has H(t) = (2/3 7 )i _1 . The 
eigenvalues of the linearization about this point are: 

Ai = -37, 

1 



A2 > 3 ~ 4 



(97 - 14) ± J (9 7 ^ 14)2 + 24(2 _ 7)(37 _ 4) ^ _ -gj^-g 



9(2- 7 )2(3 7 -4), 

For 1 < 7 < I this point exists for all values of P and is locally a sink, whence 
4P 2 (P 2 - 1) 

Over the interval | < 7 < 2 the equilibrium point becomes a local saddle and flu > 
whenever the point exists. The stability of this point is summarized in Tabic ^ 
Bianchi I brane-world model (BIb r ): 

P(P 2 -l) P 2 (P 2 -1) 2 

s = — -, A « = o, n x = i-— — { —— — '- — -, = 37-1. 

3(7-1) 9( 7 -l) 2 (3 7 -2)' 
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This equilibrium point represents an anisotropic Bianchi I brane-world cosmology since 
Q\ ^ 0. The point is physical only when £l\ > and thus P and 7 must satisfy 

P 2 (P 2 - I) 2 

(40) 

This constraint is essentially an inverted cubic function of P 2 with up to three roots. 
It can be shown that P 2 = 1 is a local max for all values of 7, and the largest 
real root is P 2 — 3j — 2. Therefore the constraint is satisfied if P 2 is less than the 
value of the largest root, i.e, P 2 < 37 — 2. In the limit as the parameter P — > 1 
this equilibrium point becomes isotropic and represents the usual Binetruy et al. 
9 Friedmann-Robertson- Walker brane-world cosmology. The exact solution at this 
equilibrium point has H{f) = (l/37)i~ 1 . The eigenvalues of the linearization about 
this equilibrium point are: 

Ai = 37, 

1 



A 2 ,3 - 2 



(97 - 7) ± J (9 7 - I) 2 - 24(7 - 1)07 - 2) (l - 2) ) 



This point is a source on the interval 1 < 7 < 2 whenever it exists. We also note the 
value 

*-- mhmh st - (41) 

The stability of this point is summarized in Tabled 
Kasner-like bulk dominated brane-world model (K^): 

£ = -P, fi = 0, ft A = 0, q=l + P 2 . 

This equilibrium point represents a Bianchi I model containing no ordinary matter 
but does contain a form of dark matter that can be interpreted as radiation as a result 
of the bulk effect U. In the limit as P — > 1 this point approaches a vacuum dominated 
Bianchi I model or Kasner model. This point exists for all values of the parameters P 
and 7. The exact solution at this equilibrium point has H(t) = (1/(2 + P 2 ))i _1 . The 
eigenvalues of the linearization about this equilibrium point are: 

Ai = 2P 2 + 4 - 37, A 2 = 2P 2 + 4 - 67, A 3 =3P 2 -1. 

This point is a source for P 2 > 3j — 2, else it is a saddle. At this equilibrium point 

n u = l-P 2 (42) 

and therefore flu < 0. The stability of this point is summarized in Tabled 
Bianchi I bulk dominated brane-world model ( e BI^): 



P + eV4 - 3P 2 



£= ^- , fi = 0, ^a = 0, q= (4-P 2 ) + eP v / 4-3P 2 . 

There are two copies of the equilibrium point, signified by the parameter e which takes 
values ±1. These equilibrium points again represent anisotropic Bianchi I models 
containing no ordinary matter but contain a form of dark radiation U as a result 
of the bulk. However, these points differ from in that + BI^ approaches the a 
general relativistic Bianchi I model containing radiation as P — > 1 while the point 
~BI^ approaches an isotropic Robertson- Walker model containing radiation. This 
point exists only if P lies in the interval 1 < P < The exact solution at this 
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Table 1. This tabic summarizes the local stability of the equilibrium points in 
the V = V = P(P 2 - 1) case. DNE=Does Not Exist 



Label 


Equilibrium Point (E, Q, Q\) 


Range of 7 


Conditions on P 


Stability 






1 ^ -4 
1 < 7 < | 


1 < P Z 


Sink 


BI P / 


( 2P(P 2 -1) -, 8P 2 (P 2 -1) 2 n \ 
\ 3(2-7) > L 9(2-7)^(3 7 -4) ' U J 


|< 7 <2 


1 < P 2 < Pi 


Saddle 






|<7<2 


Pi < P 2 


DNE 


BIb r 


( P(P--l) „ 1 P 2 (P 2 -1) 2 \ 
\ 3(7-1) > ' 9(7-l) 2 (37-2) } 


1 < 7 < 2 


P 2 < 37 - 2 


Source 






1 < 7 < 2 


P 2 > 37 - 2 


DNE 




(-P,0,0) 


1 < 7 < 2 


37 - 2 < P 2 
1 < P 2 < 37 - 2 


Source 
Saddle 


+ Bl u 


(±(P + V4-3P 2 ),0,o) 


1 < 7 < 2 


1< P 2 < | 


Saddle 






1 < 7 < 2 


| < P 2 ' 


DNE 




(±(P-V4-3P 2 ),0,o) 


1<7<| 

4 < ~ < 14 
3 ^ ' ^ 9 


1 < P 2 < f 
1 < P 2 < Pi 


Saddle 
Sink 


-Blu 


3 v ' v 9 


Pi < P 2 < | 


Saddle 




T<7<2 


1< P 2 < \ 


Sink 






1 < 7 < 2 


4 <P 2 
3 ^ 


DNE 



where P x = \ ((8 - 3 7 ) + V(20 - 9 7 )(3 7 - 4)) 



equilibrium point has H{f) = (l/(q + l))t . The eigenvalues of the linearization 
about this equilibrium point are: 

4 - 3P 2 + 3ePV4 - 3P 2 

^1 — n ' 



A 2 = 6 - P 2 + ePy/A - 3P 2 - 37, 
A 3 = 6 - P 2 + eP ^4 - 3P 2 - 67. 
Whenever this point exists this point is a saddle. The remaining variable 

o r>4 

VL U = — > 0. (43) 
The stability of these points is summarized in Tabled 

4-1-3. Observations We observe that in general the dynamics in the case of V = Vq 
a constant, are quite different than what is found currently in the literature. However, 
as we have observed in the preceding analysis, the qualitative behaviour approaches 
the qualitative behaviour discovered by Campos and Sopuerta |15lllfi| in their analysis 
of the Bianchi I brane-worlds containing a perfect fluid. In our analysis, for P 7^ 1 
(Vq 7^ 0) we have found that none of the equilibrium points are isotropic, and all have 
some form of dark matter in the form of radiation, that is hi 7^ 0. It appears that 
the non-local bulk corrections have a significant impact on the behaviour. The initial 
singularity represented by either BIf, r and Ku are not isotropic. This result suggests 
that the initial singularity in braneworld cosmologies is not isotropic as argued by 
Coley in As stated earlier, at each of the equilibrium points here, we have 

hi = Uot~ 2 and — [P^ot^ 2 where hi is a constant and [P^Jo is a matrix of 
constants. 
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Of course the ansatz that 



V 



-rr V^ v 



« 2 A 



in- nf bl ' 
2 u iiv u 



a constant may not be physically plausible, however, it is no less acceptable than the 
alternative of artificially setting = or V^ v a liv = apriori. In the next two 
subsections we consider V to be dynamic and analyze two alternative possibilities. 



12. V =V^ a , a > 

4-2.1. Equations To illustrate this case we let a = 1 and again use equation (|34J) 
to eliminate the variable VLjj from the dynamical system (!.''>( Ill whence we obtain the 
following three dimensional system 

£' =£(g-2)+7> s £, 

9! = n[2q- (37-2)], (44) 
0^ = fi A [2g-(6 7 -2)], 

where again 

q = 1 + S 2 + i(3 7 - 4)fi + i(67 - 4)f! A 

and Vt, € K- We note that S = is an invariant set and that the dynamical 
system (|44|l is invariant under the transformation S — > — E. Therefore without loss 
of generality we restrict ourselves to the invariant set S > 0. The equilibrium points 
and their local stability are discussed below. 



4-2.2. Equilibrium Points FRW perfect fluid model (F p f): 

s = o, n = i, n A = o, q = ?2_l. 

This equilibrium point represents a zero curvature isotropic Friedman-Robcrtson- 
Walker perfect fluid cosmology and in the context of brane-world cosmologies when 
extended to the bulk is represented by the solutions in ^3 ^] . The exact solution 
at this equilibrium point has H(t) = (2/37)t _1 .The eigenvalues of the linearization 
about this equilibrium point are: 

Ai = -3 7 , A 2 = 3 7 -4, A 3 =Ps + ^(7-2). 

This point is a sink if 1 < 7 < | and Vy. < — §(7 — 2), else this point is a saddle. At 
this point flu = and the stability of this point is summarized in Tabled 
FRW brane-world model (Fb r ): 

e = o, n = o, Oa = i, 9 = 37-1. 

This equilibrium point represents the Binetruy et al. 9 Friedmann-Robcrtson- 
Walker brane-world solution. The exact solution at this equilibrium point has 
H(t) — (l/37)t~ 1 . The eigenvalues of the linearization about this equilibrium point 
are: 



Ai=3 7 , A 2 =6 7 -4, A 3 =Ps+ 3(7-1). 
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Table 2. This tabic summarizes the local stability of the equilibrium points in 
the V = V S T, case. DNE=Does Not Exist 



Label 


Equilibrium Point (S, 0, 


Range of 7 


Conditions on Vt, 


Stability 


F pf 


(0,1,0) 


1< 7 < | 
1 < 7 < I 
l<7<2 


p E <-|( 7 -2) 
p E >-|( 7 -2) 


Sink 
Saddle 
Saddle 


Fbr 


(0,0,1) 


1 < 7 < 2 


< -3(7 - 1) 
> -3(7 - 1) 


Saddle 
Source 


F u 


(0,0,0) 


1<7<| 
I<7<2 


Pe e k 
Pe < 1 


Saddle 
Sink 




(Vi-P s ,o,o) 


1 < 7 < 2 


Pe < -3(7 - 1) 
-3(7-1) 
Pn > 1 


Source 
Saddle 
DNE 



This point is a source for 1 < 7 < 2 and "Ps > —3(7 — 1), else this point is a saddle. 
At this point flu — and the stability of this point is summarized in Table [3 
FRW bulk dominated brane-world model (F^): 

s = 0, = 0, n x = 

This point represents a spatially homogeneous and isotropic cosmological model 
containing no ordinary matter but does contain a form of dark matter in the form 
of radiation as a result of the bulk term U. The eigenvalues of the linearization about 
this equilibrium point are: 

Ai=Pe-1, A 2 =4-3 7 , A 3 =4-67, 

where 

n u = 1. 

Over the interval | < 7 < 2 when Vs < 1 this point is a sink, else the point is a 
saddle. The stability of this point is summarized in Table [3 
Kasner-like bulk dominated brane-world model (K^): 

n = y/i-Vz, n = o, r>A = o, q = 2-v^. 

This equilibrium point represents a Bianchi I model containing no ordinary matter 
but does contain a form dark matter in the form of radiation as a result of the bulk 
term lA. This point is physical only when < 1. In the limit as ^ this point 
approaches a vacuum dominated Bianchi I model or Kasner model. The eigenvalues 
of the linearization about this equilibrium point are: 

Ai=2-2-Ps, A 2 = 6 - 2-Pe - 3 7 , A 3 = 6 - 2-p s - 67 

where 

This point is a source whenever T^s < —3(7 — 1), else it is a saddle. The stability of 
this point is summarized in Table [21 
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4-2.3. Observations The most interesting result in this case is that the Binetruy et 
al. [HI Friedmann-Robertson- Walker brane-world solution may not represent the past 
asymptotic state. In particular there are parameter values of such that the initial 
singularity is represented by the somewhat more conventional Kasner-like solution. We 
also find these cosmologies do not necessarily evolve to a purely general relativistic 
cosmology since one of the local sinks is dominated by the bulk term U. 

4-3- V = VnuS£,P>0 

4-3.1. Equations Similar to the previous case we will set (3=1 for illustrative 
purposes only and again use equation l|34[l to eliminate the variable flu from the 
dynamical system (|30ll . We thus obtain the following three dimensional system 

s' =£(<?- 2) + nai-£ 2 -n-cix), 

Q' = Q[2q- (37-2)], (45) 

n' A = n A [2 ? -(67-2)], 

where again 

q = 1 + S 2 + ~(3 7 - 4)n + ~(6 7 - 4)fi A , 
and Vn u £ M. The equilibrium points and their local stability are discussed below. 



4-3.2. Equilibrium Points FRW perfect fluid model (F p f): 

e = o, Q = i, n x = o, q = ^Y^ 

This equilibrium point exists for all parameter values and represents a zero curvature 
isotropic Friedman-Robertson- Walker perfect fluid cosmology and in the context of 
brane-world cosmologies when extended to the bulk is represented by the solutions 
in |2B[n]. The exact solution at this equilibrium point has H{t) = (2/37)< _1 . The 
eigenvalues for this point are: 

3 

Ai = -37, A 2 = 37-4, A 3 = - 7 -3. 

This point is a sink over the interval 1 < 7 < |, else it is a saddle. At this point 
Vlu = 0. The stability of this point is summarized in Table |3| 
FRW brane-world model (F^): 

E = 0, O = 0, A = 1, g = 3 7 -l. 

This equilibrium point also exists for all parameter values and represents the the 
Binetruy-Deffayet-Langlois zero curvature brane-world solution ■ The exact solution 
at this equilibrium point has H(f) = (l/37)i _1 . The eigenvalues for this point are: 

Ai = 37, A 2 = 67 - 4, A 3 = 37 - 3. 

This point is a source over the interval 1 < 7 < 2, where again ilu = 0. The stability 
of this point is summarized in Table [31 
Bianchi I vacuum model (Kasner) ( e K): 



E = e, O = 0, ft A = 0, q = 2. 
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The parameter e = ±1 signifying that there are two of these points that exist for 
all values of the parameters 7,'Pn M . These equilibrium points represent anisotropic 
vacuum Bianchi I models, known as Kasner cosmological models. The exact solution 
at this equilibrium point has H(t) = (l/3)i -1 . The eigenvalues for these point are: 

X 1 =2-2eVn^ A 2 = 6-3 7 , A 3 = 6 - 67. 

This point is a saddle over the interval 1 < 7 < 2, regardless of the value of Vu u - At 
this equilibrium point Qu — 0. The stability of this point is summarized in Table |3J 
Kasner-like bulk dominated brane-world model (K^): 

E = V Qu , O = 0, O A = 0, q=l + Vl u . 

This equilibrium point represents a Bianchi I model containing no ordinary matter 
but does contain a form of dark radiation Wasa result of the bulk. The exact solution 
at this equilibrium point has H(f) = (l/'Po M )^ 1 - The eigenvalues for this point are: 

Ai=7>n w -1, A 2 = 27^+4-37, A 3 = 2p£ u +4-6 7 . 

This point is a source over the interval 1 < 7 < 2, given that Vf 2u > 37 - 2. On the 
interval | < 7 < 2 and given that < §7 — 2 this point behaves as a sink. For any 
other parameter values this point is a saddle. At this point 

Vu = l-Vl u . (46) 

The stability of this point is summarized in Table |31 
Bianchi I perfect fluid model (BI p /): 



3 7 -4 _ 1 3 7 ^4 _ n _3 7 -2 



fl = l ^— , Oa = 0, q 



This equilibrium point represents an anisotropic Bianchi I perfect fluid cosmology 
This point is physical when Q > 0, which implies that Vn u and 7 must satisfy 
3 7 -4 

The eigenvalues for this point are: 
Ai = -37, 



1 - ~^r- > 0- (47) 
2V? 



where 



*2. 3 = i U( 7 - 2) ± J[3( 7 - i)f + 24(3 T - - 2) (l - 

n „ _ ^-4)(2- 7 ) (4g) 
This point, whenever it exists is a saddle over the interval 1 < 7 < |, and it is a sink 



4 
3 

Bianchi I brane-world model (BIf, r ) 



on the interval | < 7 < 2. The stability of this point is summarized in Table |3J 



37-2 37-2 

1 &U 1 flu 
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Table 3. This tabic summarizes the local stability of the equilibrium points in 
the V = Vn u &u case. DNE=Does Not Exist 



Label 


Equilibrium Point (£, 0, ilx) 


Range of 7 


Conditions on T 5 ^ 


Stability 


b pf 


(U, 1, UJ 


1 < 7 < 3 

3 < 7 < z 


V 2 p. K+ 


Sink 
Saddlp 


Fbr 


(0,0,1) 


1 < 7 < 2 




Source 


e K 


(e,0,0) 


1 < 7 < 2 


T'L G » + 


Saddle 






1 < 7 < 2 




Source 


K u 


(^,0,0) 


1<7<| 
| <7<2 


ni<3 7 -2 

§7-2<PL<37-2 

nl, < ?7 - 2 


Saddle 
Saddle 






|<7<2 


Sink 






1<7<| 


^ G M+ 


Saddle 


Bl pf 


/ 3 7 -4 , 3 7 -4 n \ 
V % ' ' J 


!<7<2 




Sink 






|< 7 <2 


^<§7-2 


DNE 


BIb r 




1 < 7 < 2 


^„>37-2 


Saddle 






1 < 7 < 2 


^„<3 7 -2 


DNE 



This equilibrium point represents an anisotropic Bianchi I brane-world cosmology since 
Q\ =/= 0. This point is physical only when £l\ > and thus Vq u and 7 must satisfy 

3 -" 2 ^ (49) 



1 



V 2 
1 Qu 



> 0. 



The eigenvalues for this point are: 
Ai = 37, 

1 

2 



A 



2,3 



3( 7 - 1) ± J [3(7 - I)] 2 + 24(3 7 - 2)( 7 - 1) 1 - 



where 



3(3 7 -2)( 7 -l) 



V 2 
1 Qu 



(37 - 2) 

V 2 
1 n u 



(50) 



This point, whenever it exists, is a saddle over the interval 1 < 7 < 2, its stability is 
independent of the parameter Vn u ■ The stability of this point is summarized in Table 

El 



4-3.3. Observations An interesting observation in this case is that the Binetruy- 
Deffayet-Langlois Braneworld solution is always a source. However, there will exist 
cosmological models that will asymptote to the past towards the anisotropic bulk- 
dominated Kasner-like models represented by the points Ku when Va u > 37 — 2. 
This means that the initial singularity is again not necessarily isotropic. The future 
behaviour is determined primarily by the parameter 7. If 7 < 4/3 then the general 
relativistic flat Robertson- Walker perfect fluid model is the future attractor. However, 
if 7 > 4/3 then the non-local bulk corrections become dominant, that is, there exists 
models in which fly -/+ as r — > 00. Depending on the parameter Vq u either K-u or 
BIpj becomes the future asymptotic state. It is worth noting that neither of these 
two local sinks are isotropic. 
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5. Conclusions 

A dimensionless variable V was introduced to model the non-local bulk corrections 
arising from the bulk term V^v That is, we introduced the quantity 

V = ^_^ 

At 2 A aH 2 1 ' 

to model the non-local bulk corrections. Since there are no equations in the brane- 
world approach to model the effects of the non-local bulk term "P M „ , there is no equation 
to determine the effects of V . Therefore, one is required to make justifiable and 
reasonable assumptions on the form of V . Here, we assume that V is a differentiable 
function of the remaining dynamical variables which results in self-similar asymptotic 
solutions. 

In this paper we specifically investigated three different cases, V = Po, V = 'PsS, 
and V = Vvtu^-u where Vq,V-e and Va u are constants. It was found that in each 
case, the effects of the non-local bulk effect were significant. In the standard 
cosmological brane-world situation, the Binetruy-Deffayet-Langlois zero curvature 
brane-world solution [!|| represents the 'generic' initial behaviour in brane-world 
cosmologies containing ordinary matter. Here we find that Kasner-like bulk-dominated 
models can also act as sources in all three cases and that the Binetruy-Deffayet- 
Langlois zero curvature brane-world solution does not even exist as an equilibrium 
state in one of the cases studied. With regards to possible future behaviors, we find 
for 1 < 7 < 4/3 that the future asymptotic state is gcncrically matter dominated, in 
that the future asymptotic state has tt p f -/-> as r — > oo. For 7 > 4/3 we find that 
there exist future asymptotic states that can become 'bulk-dominated' in that 
as r — > 00. 

The qualitative analysis of this system is by no means complete. Only a local 
analysis of finite equilibrium points was investigated, none of the global behaviour 
was determined and none of the bifurcations were analyzed in any detail. There 
may be some additional interesting dynamics to be found, at points of infinity or 
even closed/periodic orbits at finite values. Further work in this area should include a 
discussion of these topics, but the motivation for this paper was to illustrate the effects 
of the non-zero non-local bulk effect V^ v in which one makes the assumption that all 
equilibrium states of the resulting dynamical system represent self-similar cosmologies. 

In conclusion, it should be apparent that neglecting the bulk effects in the brane- 
world approach to brane-world cosmologies is extremely dangerous. Even in this 
simple toy model, with a relatively simple ansatz for the non-local bulk term V^, 
we find significant differences with what is currently found in the literature when one 
ignores these effects. 
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